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Abstract 

The solutions of a renormalized BCS model are studied in two space dimensions in 
s, p and d waves for finite-range separable potentials. The gap parameter, the critical 
temperature T c , the coherence length £ and the jump in specific heat at T c as a function 
of zero-temperature condensation energy exhibit universal scalings. In the weak-coupling 
limit, the present model yields a small £ and large T c appropriate to those for high-T c 
cuprates. The specific heat, penetration depth and thermal conductivity as a function of 
temperature show universal scaling in p and d waves. 

PACS number(s): 74.20.Fg, 74.72.-h, 74.76.-w 

1 Introduction 

At low temperature and in the weak coupling limit, a collection of weakly interacting electrons, 
spontaneously form large overlapping Cooper pairs |l| leading to the microscopic Bardeen- 
Cooper-Schreiffer (BCS) theory of superconductivity |2|, [3]. For usual superconductors, the 
s-wave BCS theory yields £kp ~ 1000 in agreement with experiments, where £ is the coherence 
length and kp the Fermi momentum. There has been renewed interest in this problem with the 
discovery of high-T c superconductors. The high-T c materials have a small £: £kp ~ 10 ||, |5fl. 
Inspite of much effort, the normal state of the high-T c superconductors has not been satisfac- 
torily understood ||. There are controversies about the appropriate microscopic hamiltonian, 
pairing mechanism, and gap parameter ||, [7J. 

Many high-T c superconductors have a conducting structure similar to a two-dimensional 
layer of carriers f|, 0, U, which suggests the use of two-dimensional models. Moreover, there 
are evidences that the high-T c cuprates have singlet <i-wave Cooper pairs and the gap parameter 
has d x 2_ y 2 symmetry in two dimensions j7|. Recent measurements of penetration depth A(T) 



and superconducting specific heat at different temperatures T [10] and related theoretical 



analyses O, [L2l also support this point of view. According to the isotropic s-wave BCS theory, 
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as T — > 0, both observables exhibit exponential dependence on T 0, |12| . The experimental 
power-law dependence of these observables on T can be explained by considering anisotropic 



gap parameter with node(s) on the Fermi surface in higher partial wave [[H], 12]. It seems that 



some properties of the high-T c cuprates can be explained by the <i-wave BCS equation in two 
dimensions for weak coupling. Yet there has been no systematic study of the BCS equation in 
non-s waves in two dimensions. 

The BCS theory considers N electrons of spacing L, interacting via a weak potential of 
short range such that r$ « L and « R where R is the pair radius. When suitably 
scaled, most properties of the system should be insensitive to the details of the potential and 
be universal functions of the dimensionless variable L/R [Tj|. The usual BCS treatment |J 
employs a phonon-induced two-electron potential of moderate range. In this paper we study 
the weak-coupling BCS problem in two dimensions for s, p, and d waves with two objectives in 
mind. The first objective is to identify the universal nature of the solution and its relation to 
to high-T c superconductors. The second objective is to find out to what extent the universal 
nature of the solution is modified in the presence of realistic finite-range (nonlocal separable) 
potentials. Instead of solving the BCS equation on the lattice with appropriate symmetry, we 
solved the equations in the continuum. This procedure should suffice for present objectives. 

In place of the phonon-induced BCS model we employ a renormalized BCS model with sep- 
arable potential which has certain advantage. The renormalized BCS model leads to convergent 
result even in the absence of potential form factors or momentum/energy cut off, as required 
in the standard BCS model. The original BCS model yields a linear correlation between T c 
and Td, where To is the Debye temperature. This correlation was fundamental in explaining 
the observed isotope effect in conventional superconductors. The high-T c materials exhibit an 
anomalously negligible isotope effect and a linear correlation between T c and Tp, where Tp 
is the Fermi temperature. This suggests a different interaction for superconductivity in the 
high-T c materials. The present renormalized model yields a linear scaling between T c and Tp. 
Because of this scaling with Tp, unlike in the phonon-induced BCS model, the present T c can 
be large and appropriate for the high-T c cuprates in the weak coupling region. The present 
model also produces an appropriate T c /Tp ratio and a small £fcp in the weak coupling region 
in accord with recent experiments || on high-T c cuprates. In spite of these we are aware that 
there are controversies in the description of the high-T c materials, for example, the microscopic 
hamiltonian and the pairing mechanism. Also, the normal state away from T c seems to be very 
different from a standard Fermi liquid. However, we find that there are certain characteristics 
of these high-T c materials which can be studied within the present renormalized mean-field 
BCS model based on the standard Fermi liquid theory. 

Previously, there have been studies of this problem in two dimensions in terms of two-body 
binding in vacuum and Cooper pair binding employing short- and zero-range potentials 0, [5], [L4| . 
Such studies have not fully revealed the universal nature of the solution. In the present paper 
we employ the zero-temperature condensation energy per particle, AU, of the BCS condensate 
as the reference variable for studying the BCS problem. As the condensation energy increases, 
one passes from the weak to medium coupling. Also, AU is a physical observable and is the 
appropriate reference variable as we shall see. In this paper we calculate the zero-temperature 
gap parameter A(0), the critical temperature T c and the specific heat per particle C S (T C ) in 
all partial waves (denoted by angular momentum I) and the zero-temperature pair size in s 
wave. We find that these observables obey robust universal scaling as functions of AU valid 



2 



over several decades in the weak to medium coupling region independent of the potential range 
parameter. 

Here, we also calculate the temperature dependence of different quantities, such as, the 
gap parameter A(T), C S (T), A(T) and thermal conductivity K{T) for T < T c . Of these, the 
T dependence of C S (T), A(T) and K(T) are specially interesting. For isotropic s wave, the 
BCS theory yields exponential dependence on temperature as T — > for these observables 
independent of space dimension || |T2]| . The observed power-law dependence of some these 
observables can be explained with anisotropic gap parameter in non-s waves with node on 
the Fermi surface. We find universal power-law dependence in both cases for non-s waves 
independent of the range of potential. For I ^ we find C S {T) w 2C n (T c )(T /T c ) 2 and K S {T) w 
K n (T)(T '/T^ 1 - 2 valid for almost the entire range of temperature. The suffix n and s refer to 
normal and superconducting states, respectively. Similar dependence was predicted from an 



analysis of experimental data |10| as well as from a calculation based on Eliashberg equation 
In order to detect the anisotropy in A(T) it is appropriate to consider the function 
AA = (A(T) - A(0))/A(0) @. We find that for small T it behaves as AA ~ (T/T c ) 13 . Similar 
power-law dependence was conjectured before . 

From the weak-coupling BCS equation we establish the following relations analytically: 
2A(0)/T C ^3.528 (3.026), A(0) = 2^A~U (2y/AU), T c « 1.13AVaU (1.322 v / AC7), AC/C n {T c ) « 
1.43 (1.05), C S {T C ) « Vtitty/W (8.915VAT7), AU/U n (T c ) « 0.473 (0.348), for I = (^ 0). 
Here AC is the jump in specific heat at T = T c , C (U) is the specific heat (internal energy) 
per particle. 

The plan of the paper is as follows. In section 2 we present the renormalized BCS model. 
In section 3 we present an account of our numerical study. Finally, in section 4 we present a 
brief summary of our findings. 



2 Renormalized BCS model 

The two-body problem, Cooper and BCS models all exhibit ultraviolet divergences for zero 
range potential and require regularization or renormalization to produce finite results PJ. The 
usual BCS model employs a cut-off regularization (Debye cut off) for obtaining finite observ- 
ables. In momentum space, the standard BCS model interaction is taken to be constant for 
momenta between 2m(Ep — Ed)/?i 2 and 2m(Ep + Ed)/?i 2 and zero elsewhere with Ed (Ef) 
the Debye (Fermi) energy. This implies a moderate range of the interaction. This potential is a 
physically motivated phonon induced electron-electron one 0. The present renormalized BCS 
equation may lead to well-defined solution without requiring a cut off even in the absence of 
potential form factors. The present model with finite-range potential also leads to a well-defined 
mathematical problem. 

We consider a two-body system, each of mass m, in the center of mass frame Q. The single 
(two) particle energy is given by e q = h 2 q 2 /2m (2e 9 ), where q is the wave number. We consider 
a purely attractive short-range separable potential in partial wave I: 

Vp q = -V cig p ig q i cos(/#) (1) 

where 9 is the angle between vectors p and q and q = 1 (2) for I = 0) as in Ref. [[0J. 
Here g pi and g q i are potential form factors and V is the potential strength. The potential V m 
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is the effective electron-electron potential in the superconductor in the presence of lattice and 
other electrons. The Schrodinger equation in this case leads to the following condition for a 
two-particle bound state in vacuum of binding B 2 is 



c^4cos 2 (W)(B 2 



^r 1 



(2) 



where 6 is the angle of vector q and e q = Ti q /2m with m the mass. In order to simplify the 
notation, the trivial I dependence of many quantities, such as Vp q , B 2 , etc., is suppressed. 

For attractive potential (|2|), two electrons on the top of the full Fermi sea at zero temperature 
shows pairing instability. Singlet (triplet) pairing is assumed to take place in even (odd) partial 
waves. The Cooper pair problem for two electrons above the filled Fermi sea for this potential 
is given by 0, § Vq 1 = ci J2q>i 9 q i cos 2 (W)(2e q - 2E)-\ with Cooper binding B c = 2 - 2E. 
Using (E3), the Cooper problem is written as 



J29 2 g iCos 2 (W)(B 2 + 2e q 



J29 2 ql cos 2 (W)(2e q -2EY 
0. 



At a finite T, one has the following BCS and number equations 

A n = 



VVpq^r-tanh^ 
V 2E * 



N 



E 



i - 



tanh 



2T' 
E„ 



2T 



(3) 

(4) 
(5) 



with Eq = [(e q — fx) 2 + [Aql 2 ] 1 / 2 , where A q is the gap function. Unless the units of the variables 
are explicitly mentioned, in (|2|) — (H) and in the following all energy/momentum variables are 
expressed in units of Ep, such that fi = fi/Ep, T = T/Tp, q = q/kp, E^ = E^/Ep, etc, where 
ji is the chemical potential. Here A q has the following anisotropic form: A q = g q iA y/c l cos(W) 
where A and g q i are dimensionless. The usual BCS gap is defined by A(T) = fl , g (=i)zA , which 
is the root-mean-square average of A q on the Fermi surface. Using the above form of A q , the 
BCS equation (01) can be written as 



1 

Vn 



9ql COS 2 '(W) 



1 



2E a 



tanh — ^. 
2T 



(6) 



Now (0) and 



lead to 



q 

The summation is evaluated according to 

E 



2e + B 2 



■ tanh q 



47T 



Er, 



N 



2T 



0. 



(7) 



de q d6 = — I de„ 
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where N is the number of electrons. As in the standard BCS model, we have a constant density 
of state in (§) but now with the generalization to include the angular dependence. With the 
help of (|8]), (|5|) and (0) can be explicitly written as 



de q d6 



1 



^ tanh — - 
2T 



An 



(9) 



de q d9g ql cos 2 (W) 



2e„ + B 2 E Q 



1 tanh^ 



2T 







(10) 



respectively. In terms of Cooper-pair binding Eq. (|T0| ) can be rewritten as 

r2n 



d6 cos 2 (W) 



de„ 



9ql 



E 



Q 2 E 
<ie„— tanh q 



2T 







Even in the absence of potential form factors g p i = 1, @, (|T0|) and ([TTJ) are well defined 
without any energy /momentum cut-off, though each part of the integral in these equations 
diverges separately. This is why this model is termed renormalized. Potential ([I]) with g p i = 1 
is the zero-range delta-function potential. The standard BCS model uses essentially the above 
potential with energy/momentum cut-off for obtaining convergence. In renormalized equation 
( |T(jl ) the usual energy/momentum cut-off and the potential strength Vq have been eliminated 
in favor of the two-body binding B 2 . Recently, the role of renormalization in nonrelativistic 
quantum mechanics has been discussed [[U| . 

Equation ( |TTD has following analytic solutions for s-wave zero-range potential (g q o = 1) 
in the weak coupling limit (/i = 1). At T = 0, A(0) = V2B~ C . At T = T c (A = 0), we 
have: T c = exp(j)y / 2B c /n pa 0.8a/Bc where 7 = 0.57722... The standard BCS model yields 



in this case T c /Tc = exp(7) v2£QVvTd pa 0.8^/ B c /Td where To is the Debye temperature. 
To illustrate the advantage of the renormalized model let us consider a specific example with 
T D = 300 K, Tp = 3000 K and B c = 10 K. We take B c < 10 K as defining the weak-coupling 
region. With B c = 10 K, the standard BCS model yields T c = 44 K, whereas the present 
renormalized model yields T c = 138 K. Hence for same coupling, the renormalized model leads 
to an enhanced T c . In the standard BCS model one has to have a much larger B c , clearly 
outside the weak-coupling region, in order to have T c > 100K. 

The universal ratio 2A(0)/T C = 2tt/ exp(7) pa 3.528 remains unchanged for s wave in three 
dimensions as well as for the trivial case of anisotropic pairing (~ exp in two dimensions 
for I 7^ [14[]. In Ref. |14j we essentially changed the potential form-factors without introducing 
any explicit angular dependence in the BCS and number equations and find that the universal 
nature of the solution was unchanged with such change. Hence we expect to extract certain 
universal properties of ( |T0D analytically for I 7^ by employing the correct angular distribution 
and no potential form factors (g q i = 1). 

Next we study ([II]) analytically for weak coupling (fi = 1, A << 1) and for g q[ = 1. At 
T = 0, (0) can be integrated to yield 



2 IT 



d9 cos 2 (W) In 



[1 + A 2 (0)q cos 2 ^)] 1 / 2 - 1 



Bo 



0. 
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which for A(0) << 1 and I ^ reduces to 



In 



QA 2 (0) 



1 r 2n 



2Bo 



IT JO 



d9cos 2 (W) lncos 2 (Zfl) = 0.3863, 



or, A(0) ~ 1.213V52- At T = T c , we again have T c = exp(7)v / 2 _ B^/7r, so that we have the 
following new universal constant for I ^ 0: 2A(0)/T C rs 3.026. 
The condensation energy per particle at T = is defined by 



AU = \U, - U n 



E 2C, 



where C 



A q 

/x). This can be evaluated straightforwardly to lead to 

1 ^ 

3Wo 



AU 



d A 2 (0) cos 2 (W)d9 



which yields A(0) = 2V AU for all Z. Using the universal relation between A(0) and T c given 
above one has T c « 1.134 V / AZ7 (1.322 VST) for / = (/ ^ 0). For all Z, U n (T c ) = tt 2 T 2 /6, so 
that AU/U n {T c ) «0.473 (0.348) for Z = & 0). 

The superconducting specific heat per particle is given by 



a 



NT 2 



EA(i-/ q ) [K-\ T 



l„rfA| x 
dT 



(12) 



where / q = 1/(1 + exp(£' q /T)). The normal specific heat C n is given by (|12D with A q = 0. The 
jump in specific heat per particle at T = T c (A(T C ) = 0), AC = [C s — C„]t c is given by || 



AC 



iVT, 



E 



q l 



/A q ^ 

/q(l-/q)^T 



(13) 



In the special case g q i = 1, the radial integral in (|T^) can be evaluated as in Ref. and we get 

dA 2 (T)" 



AC = - q 



de q d6 
4ttT, 



/q(l " /q) 



cos 2 (Z#). 



(14) 



This leads to AC = — (l/2)[dA 2 (T)/dT]T=T c for all /. In a systematic (numerical) study 
we find A(T)/A(0) = B(l - T '/T c ) 1 ' 2 valid for T « T c , with 5 « 1.74 (1.70) for / = (^ 0). 
Using the value of B and the universal ratio 2A(0)/T C , we obtain AC/C n (T c ) « 1.43 (1.00) 
for Z = (^ 0), where C n (T) = vr 2 T/3. Consequently, C B {T c )/y/~MJ « 9.065 (8.915) and 
T c C n (T c )/AU « 4.229 (5.749) for Z = (^ 0). 
The penetration depth A is defined by [|J 



A~ 2 (T) = A" 2 (0) 



^ E /q(^ /q) 



The thermal conductivity ratio K s (T)/K n (T) is defined by [[17] 

K S (T) _ S q C^q/q(l-/q) 
A»CT) E q C 2 /q(l-/q) ' 



(15) 



(16) 



The denominator of ( |I6|) essentially corresponds to the normal-state thermal conductivity with 
the BCS gap set equal to zero. 

The dimensionless s-wave pair radius defined by £ 2 = (ip q \r 2 \ip q ) / {ip q \ip q ) with the pair wave 
function ip q = g q iA/ (2Eq) can be evaluated by using r 2 = —V 2 . In the weak coupling limit, 
the zero-range analytic result of Ref. § leads to £ 2 = 0.5A~ 2 (0) = 0.125(Af/)- 1 . 

3 Numerical Results 




AU FIG1 

Fig. 1 The plot of specific heat C S (T C ) (dashed line), T c (dotted line), gap parameter A(0) (dashed- 
dotted line) for s, p, and d waves and s-wave pair radius £ 2 (solid line) versus zero-temperature 
condensation energy per particle AU for different potential parameters between a = 1 to oo. For the 
first two variables there are two distinct lines; the upper one for p and d waves and the lower one for 
s wave. 

We solve coupled equations (§) and ( |10]) numerically using dimensionless form factors g q i = 
(e q ) l ^ 2 [a/(e q + a)]^ +1 ^ 2 with the correct threshold behavior for small momenta, where a is the 
range parameter. Following Ref. [§] we calculate the dimensionless gap parameter A(0) = 
5f 9 ( =1 )^A, T c , C S (T C ), the s-wave pair radius £ 2 at T = as well as A(T), A(T), and C(T) for 
different coupling. In figure 1 we plot A(0), T c , C S (T C ), and £ 2 at T = versus AU and find 
universal scalings. The calculations were repeated for different potential ranges a. We varied 
a from 1 to oo and found figure 1 to be insensitive to this variation in each partial wave. For p 
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and d waves @ and (pTO|) diverge for a — >■ oo and calculations were performed for a = 1 to 10. 
The increase in AU of figure 1 corresponds to an increase in coupling. We could express this 
change in coupling by a change in B 2 or in Cooper pair binding and plot the variables of figure 
1 in terms of these bindings as in Ref. |L4"| |. Then each value of range parameter leads to a 
distinct curve. However, if we express the variation in coupling by a variation of an observable 
of the superconductor, such as AU or T c , universal potential-independent scalings are obtained. 
In each case the exponent and the prefactor of each scaling relation are in excellent agreement 
with the analytic relation obtained above without form factors. 

2.5 
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O 
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0.5 - 



0.0 




0.0 



c FIG 2 

Fig. 2 Specific heat C s (T)/C n (T c ) versus T/T c for s (dashed line), p and d (solid line) waves and 
potential parameters between a = 1 to oo. The analytic fit (dashed-dotted line) C s (T)/C n (T c ) = 
2(T/T C ) 2 , top and d waves is also shown. 

The T c should not arbitrarily increase with coupling as figure 1 may imply. With increased 
coupling the electron pairs should form composite bosons which may undergo a phase transition 
under the action of a residual interaction. According to a numerical study this transition 
happens at a temperature of 0.1 



This is why the T c curve in figure 1 has been plotted till 
about T c = 0.1. For a very large class of two-dimensional high-T c materials, the T c has been 
estimated to be about 0.05 ||, which corresponds, for g qi = 1, to B 2 = 0.004. This small value 
of B2 is in the weak coupling region where the universality of the present study should hold. 
The corresponding dimensionless pair radius (£ 2 ~ 80) at T c = 0.05, as obtained from figure 1, 
is in agreement with experimental analysis ||. Hence, the present study is relevant for these 
high-T c materials. 
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Next we studied the temperature dependence of A(T), C S (T), A(T) and K(T) for T < T c . 
For BCS superconductors, these observables have exponential dependence on T as T — » [[3j |i~2|l . 
The two-dimensional high-T c superconductors have a power-law dependence on T. In figures 2, 
3, and 4 we plot C 8 (T)/C n (T e ), AA(T) = (A(T) - A(0))/A(0), and K s {T)/K n {T) versus T/T c , 
respectively. In these figures we find universal power-law dependence, essentially independent 
of potential range, for 1^0. We find C S {T) « 2C n (T c )(T/T c ) 2 , K S {T) « 2K n {T)(T /T c ) 1 - 2 
for almost the entire temperature range and AA(T)/A(0) ~ (T/T C ) L3 for small T/T c . The 
T 2 dependence of C S (T) was found in a theoretical study of the Eliashberg equation |yj and 
in an analysis of experimental data []10[ . The power-law dependence of A(T) on T was also 
conjectured before ||12|| . The gap function A(T) has essentially the same universal behavior as 
in s wave [|J. 
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Fig. 3 Penetration depth AA(T) versus T/T c for s (dashed line), p and d (solid line) waves and 
potential parameters between a = 1 to oo. 



4 Summary 

Scalings are established among A(0), T c , C S (T C ), and £ 2 , as function of AU, independent of 
potential range in s, p and d waves from a study of a renormalized BCS equation in two 
dimensions. The present renormalized model yields a large T c and a small £ appropriate for 
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some high-T c superconductors. The T dependence of AA(T), C S (T) and K S (T) below T c in non-s 
waves show power-law scalings distinct to some high-T c materials at low energies. No power-law 
T dependence is found in s wave for AA(T) and C S (T). Calculations performed with cos(Z6>) 
and sin(/6 ) ) angular dependences yielded identical results. Though we have used a separable 
potential, in view of the universal nature of the study we do not believe the present conclusions 
to be so peculiar as to have no general validity. We have repeated the s wave calculations 
with local Yukawa potential and found the results to be independent of potential in the weak- 
coupling region. This is also in agreement with a suggestion by Leggett [Tj|. Although, there 
are controversies about a microscopic formulation of high-T c superconductors, it seems that 
the two-dimensional (i-wave BCS equation for weak coupling can be used to explain some of 
their universal scalings. A similar study of universality has recently been performed in three 
dimensions employing the renormalized BCS equation |19[ , which also leads to an enhanced T c 
in the weak-coupling limit. 




Fig. 4 Thermal conductivity ratio K S {T) / K n {T) versus T/T c for s (dashed line), p and d (solid 
line) waves and potential parameters between a = 1 to do. 
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